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Abstract 

By means of a Madelung decomposition, exact periodic traveling so- 
lutions are constructed for a modified nonlinear Schrodinger equation 
derived by Stenflo and Gradov, describing electrostatic surface waves 
in semi-infinite plasma. The condition for the existence of bi-stable 
equilibria is discussed. A conservation law as well as the modulational 
instability admitted by the model are analyzed. 
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I Introduction 

The description of nonlinear plasma surface waves is a growing subfield in 
plasma physics [TJ [2] . Reasons for this research can be found in industrial ap- 
plications [3], plasma microwave diagnostics [I], some aspects in laser fusion 
[5] and space plasma physics [6]. In particular, much attention has been paid 
to electrostatic soliton-like oscillations, or oscillons for cylindrical [H] 

as well as spherical [101 [II] plasmas. Similar oscillon solutions have been also 
reported as highly localized ring structures in Bose-Einstein condensates [T2] . 
The present work, on the other hand, consider one- dimensional nonlinear pe- 
riodic structures in a plasma interface. These oscillatory solutions are found 
through a different treatment of a modified nonlinear Schrodinger equation 
derived earlier by Stenflo and Gradov [131 IH] • This equation resembles some 
classes of PDEs discussed in the context of dissipative quantum mechanics 
[15]- [17] and can be shown to admit exact solitary solutions [HI [T9] . 

This work is organized as follows. In Section II the generalized nonlin- 
ear Schrodinger equation for surface waves in a sharp plasma interface is 
reviewed. Instead of removing the time-derivative term as in previous contri- 
butions [13], here a Madelung transformation is applied, so that a system of 
partial differential equations for the amplitude and the phase of the envelope 
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electrostatic field is derived. In Section III, the problem of finding traveling 
solutions for this set of partial differential equations is cast in the form of an 
autonomous Hamiltonian system. Some basic properties of the model equa- 
tion are considered in Section IV, including the study of its modulational 
instability and the derivation of a conservation law. The shape of the asso- 
ciated Sagdeev potential and the different possible equilibria are analyzed. 
Illustrative trajectories are computed numerically. Section V is reserved to 
the conclusions. 



II Madelung transformation 

Our starting point is Eq. (5) of Ref. [T3], specialized to the case of a plane 
boundary (z — 0), 

— d t <p + d 2 x cf> - + 2/3 7 £#o| 4>o + /3 2 |0o| 2 0o + 2 k 2 A4> = . 

(1) 

Formerly [13], a similar equation was considered, but without the time-de- 
rivative term, as well as the terms proportional to 7 (defined below Eq. 
01])) and A = 2 — oj^/oj 2 , a measure of the mismatch to the linear surface 
plasmon frequency uj p J s/2, where u p is the plasma frequency. The purpose of 
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the present communication is to show that a plasma surface can admit exact 
one-dimensional oscillatory structures, provided these contributions are taken 
into account. In the former case, an exact localized, surface solitary wave 
was found [2]. 

Let us briefly examine the meaning of the different terms in Eq. ([!]), as 
well as its derivation [TBI [T4"] . Consider a plasma boundary at z = 0, with a 
fixed ionic background of constant density for z > and vacuum for z < 0, 
and electrostatic oscillations. Then the scalar potential 4>(x,z,t) exp[i(kx — 
oot)} satisfy the Laplace equation for z < 0, while inside the plasma it satisfies 
V • (eV0) = 0. Here the dielectric function e is 



InEq. (2),the~ (3 2 term is a small polynomial nonlinear correction. We con- 
sider the case of cold electrons and immobile ions, so that (3 = —Aek 3 /muj 2 , 
where — e and m are the electron charge and mass [2]. Further, assume 



for z > 0, where 0o is slowly varying both in space and time. Matching Eq. 
([3]) to the appropriated analytic continuation solution outside the plasma, 
taking into account the continuity of <p and the normal component of the 
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electric field at z = 0, one derive Eq. ([!]). For this purpose, one uses the first 
order Taylor expansion of k z <C k in powers of (f) and 



which implicitly define the phenomenological non-dimensional parameter 7. 
See Refs. [13], [H] for more details. 

In Ref. [13], instead of the plane xy the boundary of the plasma was 
taken as deformable due to the ponderomotive force. The plasma interface 
was parametrized by a suitable function z = £(x, t), allowing to eliminate the 
time-derivative term in Eq. ([I]) in a moving frame. Here, we follow a different 
approach, taking a fixed z = plasma interface and observing that Eq. ([I]) 
is a modified nonlinear Schrodinger equation. This suggest the application 
of a Madelung, or eikonal, decomposition, 



where the amplitude A = A(x, t) and the phase S = S(x, t) are real functions. 
Inserting Eq. ([5]) into Eq. ([1]) and separating the real and imaginary parts, 
it follows that 





0o = Aexp(iS) , 



(5) 



8 k 2 



dls 



d t S 



8 k 2 d t A 
d 2 x A [d x A] 2 



+ 2/?7 kA + P 2 A 2 + 2k 2 A . 



(6) 



(7) 



A A 2 
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Eliminating S between Eqs. (EKD, ^ follows that 

d 2 x (flgV + 2/57 fc A)^ + P 2 Al e 2 ^) + (8 fc 2 /^) 2 = , (8) 
where we have defined 

ifj = \n(A/A ) , (9) 
in terms of a reference amplitude A Q . 

Ill Sagdeev potential 

Assume a traveling wave solution ip = ip(x — ut) for Eq. ([8]), where u is 
a constant parameter. Using decaying boundary conditions, Eq. is con- 
verted into the autonomous one-dimensional Hamiltonian system d 2 ip / dX 2 = 
—dV/d^, where X = x — ut and V = V(ijj) is the Sagdeev potential 

V = 32^^j ip 2 + 2P 1 kA (e^ -^-l) + ^t(e 2 ^ -2^-1). (10) 

Hence evidently the problem is reducible to a quadrature using the energy in- 
tegral H = {jp') 2 /2 + V(i/}). Notice that the above Sagdeev potential is differ- 
ent from the one derived for similar nonlinear surface plasma wave problems 
[20l [2Tj . due to the different starting assumptions and modeling. 
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To complete the solution, Eqs. (jMZ]) can be used to obtain the phase 
S = S(X) through 

dS 8k 2 uib 00 .„„ , . ^9 l9 „,9»n 

Because the envelope 0o is slowly varying, necessarily u is a small quantity. 
However, we exclude the case where it is strictly zero. In this way, V ~ if) 2 
as if) — > — oo. Since V ~ e 2 ^ as ^ — > oo, we then have an asymmetric 
trap sustaining periodic nonlinear oscillations of if). Moreover, observing 
that V(0) = V'(0) = 0, we conclude that a double well is present whenever 
V"(0) < 0, see Fig. ([I]). In the present work we assume that the condition 
V"(0) < 0, or 

-/3 7 ^ >32^ 2 + f3 2 Al (12) 

is fulfilled, so that bi-stability is possible for a negative pseudo-energy H. 
The corresponding phase-space curves are in Fig. (|2J). 

Since (3 k is a negative quantity, to satisfy the inequality (|T2l a necessary 
condition is that 7 > 0. In view of Eqs. ([HJSD , this is a most natural require- 
ment, since we expect a bounded scalar potential. Moreover, approximating 
uj by its linear expression uj p / a/2 and using the definition of j3, it is possible 
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Figure 1: (Upper) Sagdeev potential in Eq. (jTUI) for 32(k 2 u/u) 2 = 
1/2, -2^k\A Q \ = 3,(3 2 A 2 /2 = 1/4; (Bottom) Sagdeev potential in Eq. ([TO]) 
for 32(k 2 u/u) 2 = -2^k\A \ = (3 2 A 2 /2 = 1/10. 
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Figure 2: Phase-space curves in a typical bi-stable equilibrium case of the 
Sagdeev potential in Eq. (HDD . Parameters: 32(k 2 u/uj) 2 = 1/2, -2f3-yk\A \ = 



3,(3 2 A 2 /2 = 1/A. 
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to rewrite Eq. (TT2|) as 



eA k 



) 



2 



7e A 



< 



o 




(13) 



mujp 



8 m 



which impose a lower bound on the wavelengths for bi-stable oscillations. In 
addition, the velocity u should not be too large. 

Although the problem was reduced to a quadrature, it can not be done 
in terms of elementary functions, due to the form of the Sagdeev potential. 
Nevertheless, the trajectories can be numerically calculated, as shown in Figs. 
to (JSJ), containing the auxiliary function ip, the amplitude A and the real 
part of the envelope 4>q. In the later, one can see the rapid oscillations due 
to the incorporation of the phase given by Eq. ( TTTT) . Figs. (J3KHJ) were build 
for a fixed bi-stable case but varying pseudo-energy H and initial position, 
so that ip oscillates around the left equilibrium point, the right equilibrium 
point or zero (for positive pseudo-energy). A significant distortion takes place 
according to the degree of nonlinearity. In particular, for positive pseudo- 
energy, the trajectories of the auxiliary function ip are not restricted to the 
left or right equilibria, and nonlinear oscillations around the origin take place, 



see Figs. <^M>- 
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Figure 3: Auxiliary function ip (on the top), amplitude A (on the mid) and 
real part of the envelope <po (on the bottom) for linear oscillations near the 
fixed point on the left for the Sagdeev potential in Eq. fTTOj) . Parameters: 
32(k 2 u/cu) 2 = l/2,-2Pjk\A \ = 3,p 2 A 2 /2 = 1/4, A = 0. Initial conditions: 
tp(0) = -2.2,7/7(0) = 0.5. 
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Figure 4: Auxiliary function ip (on the top), amplitude A (on the mid) and 

real part of the envelope 0o ( on the bottom) for nonlinear oscillations around 

the fixed point on the left for the Sagdeev potential in Eq. f[TU|) . Parameters: 

32(k 2 u/uj) 2 = 1/2, -2p-fk\A \ = 3,f3 2 A 2 Q /2 = 1/4, A = 0. Initial conditions: 

^(0) = -2.2,?/y(0) = 1.2. The pseudo-energy is then H = -10~ 3 . 
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Figure 5: Auxiliary function ip (on the top), amplitude A (on the mid) and 
real part of the envelope <fto (on the bottom) for linear oscillations near the 
fixed point on the right of the Sagdeev potential. Parameters: 32(k 2 u/uj) 2 = 
l/2,-2(3~fk\A \ = 3,(3 2 A 2 /2 = 1/4, A = 0. Initial conditions: ^(0) = 
1.4, ^'(0) = 0.5. The pseudo-energy is then H = —0.7, near the bottom of 
the right potential well. 
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Figure 6: Auxiliary function ip (on the top), amplitude A (on the mid) 

and real part of the envelope 0o (on the bottom) for nonlinear oscillations 

around the fixed point on the right of the Sagdeev potential. Parameters: 

32(k 2 u/uj) 2 = 1/2, -2p-fk\A \ = 3,(3 2 A 2 /2 = 1/4, A = 0. Initial conditions: 

■^(0) = 1.4, ^'(0) = 1.3. The pseudo-energy is then H = — 1CT 5 . 
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Figure 7: Auxiliary function ip (on the top), amplitude A (on the mid) 

and real part of the envelope 0o (on the bottom) for oscillations around 

zero, for a slightly positive pseudo-energy. Parameters: 32(k 2 u/u) 2 = 

1/2, -2/^7 A;|Ao| = 3,/3 2 A 2 /2 = 1/4, A = 0. Initial conditions: ^(0) = 

0, ^'(0) = y/2/10. The pseudo-energy is then H = 1CT 2 . 
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Figure 8: Auxiliary function ip (on the top), amplitude A (on the mid) 
and real part of the envelope O ( on the bottom) for oscillations around 
zero, for large pseudo-energy. Parameters: 32(k 2 u/uj) 2 = 1/2, —2/3'y k\Ao\ = 
3,p 2 A 2 /2 = 1/4, A = 0. Initial conditions: ^(0) = O,-0'(O) = y/2. The 
pseudo-energy is then H = 1. 
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IV A conservation law for the Gradov-Stenflo 
equation and modulational instability 

A peculiar property of Eq. ([1]) is the singular character of its third term when 
the envelope (fi approaches zero at some point, except in the non-generic 
situation where the derivative is also zero at that point. The existence of a 
similar singularity was even an argument to discard a dissipative nonlinear 
quantum mechanics model introduced by Kibble [22] . Actually, it is apparent 
from the original derivation [T3J E] that a formal expansion in powers of 
/c _1 <91n (f>o/dx was used. Clearly such a formal expansion is not valid near a 
singularity. Therefore, the whole theory becomes invalid in such a case. 

In this context, it is natural to assume an slowly varying and everywhere 
nonzero initial condition, (f>o(x, 0) ^ 0, \d x (j>o(x, 0)/4>o(x, 0)| <C k, Vx. A obvi- 
ous question is then: is the nonzero property of the envelope field preserved 
by Eq. ([T])? Fortunately, the answer is yes, due to the following argument. 

By inspection of Eq. a conservation law is easily detected, 



for e.g. decaying or periodic boundary conditions. Therefore, we have a con- 
straint preventing the singularity, provided the initial envelope is everywhere 




(14) 
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nonzero. 



Also due to the nonzero property of the envelope, it is convenient to 
linearize Eq. (fT]) around the exact homogeneous solution 



) = 0o exp 



8k 2 



(/3 2 |0or + 2/?7^l0olM 



(15) 



representing a plane wave whose frequency is nonlinearly shifted, where 0o 
is a positive amplitude. For simplicity, the case of zero frequency mismatch 
is being considered (A = 0). 

Assuming a perturbation in the form 



exp 



8k 2 



{(3 1 \M I + 2^k\<p \)uot + i6 



(16) 



where and 8 are real first-order quantities proportional to exp(ikx + at), 
we derive the dispersion relation 



<jj k 
8k 2 



2/? 7 £# o | + 2/3 2 |0 o r-£; 



(17) 



Notice that the potentially singular term of Eq. ([T]) has no role when lin- 
earizing around a nonzero solution. In this context this term represents a 
higher-order nonlinearity. 



From Eq. (TlTl) . it follows that a is imaginary and the wave amplitude 
remains bounded when |0 O | < —^kj(3 > 0, which is possible for k > since 

j//3 < 0. 
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On the other hand, a 2 > when |0 O | > -7 k/fi and 2/^7 k\4> \+2 /3 2 \4> \ 2 > 
k 2 , so that the amplitude is exponentially amplified. This is a modified 
modulational instability admitted by Eq. ([T]). 

V Conclusion 

As one of our final remarks, we observe that apparently there is no model 
for the parameter 7 in the literature, in terms of physical quantities. At 
the present stage, it should be considered as a phenomenological constant, 
introduced as a result of a Taylor expansion, see Eq. (J4]). Also the velocity 
u can be arbitrarily chosen. Other improvements could be the consideration 
of moving or diffuse j2] boundaries. As discussed in Section IV, the model 
Eq. (1) was derived on the assumption of a slowly varying nonzero envelope. 
It would be an relevant improvement, to derive the corresponding evolution 
equation allowing for zero amplitude solutions. These open questions cer- 
tainly deserve more attention. 
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